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First-order transition in Potts models with "invisible" states 

Rigorous proofs 

Aernout C. D. van Enter,*) Giulio Iacobelli,**) Siamak Taati 

Johann Bernoulli Institute, Rijksuniversiteit Groningen, Nijenborgh 9, 
974 7 AG Groningen, The Netherlands 

In some recent papers by Tamura, Tanaka and Kawashimap^^^ a class of Potts models 
with "invisible" states was introduced, for which the authors argued by numerical arguments 
and by a mean-field analysis that a first-order transition occurs. Here we show that the exis- 
tence of this first-order transition can be proven rigorously, by relatively minor adaptations 
of existing proofs for ordinary Potts models. In our argument we present a random-cluster 
representation for the model, which might be of independent interest. 

§1. Introduction 

In Refs. [20l) . l2ip . the authors introduced a class of Potts models, in which next 
to q ordinary -visible- colours (the Potts states), between which a standard ferro- 
magnetic nearest-neighbour Potts interaction exists, r "invisible" colours (states) 
are possible, which have zero interaction energy with any neighbour, whatever state 
that neighbour is in. 

Although the number of ground states and low-temperature states equals q, 
and there is at low temperatures spontaneous symmetry breaking of the g-fold per- 
mutation symmetry just as in the standard g-state Potts model, the transition for 
low q = 2,3,4 and high r is different from the second-order transition of the or- 
dinary two-dimensional g-state Potts model. In fact a first-order transition in the 
temperature-parameter appears. 

The occurrence of such a first-order transition contradicts a simple form of uni- 
versality which would predict that all systems with the same broken symmetry in 
the same dimension with short-range interactions have the same type of transition. 

However, such a universality property is known to be too strong to be true. 
The question of first-order versus second-order is not a universal question. Some 
counterexamples illustrating this point are the two-dimensional 3-state Kac-Potts 
model, in which a first-order transition occurs in presence of a broken 3- fold rota- 
tion (= permutation) symmetry, or the three-dimensional versions of the nonlinear 
0(n)-models treated in Ref . [1]) ,2]) ,[7])-[9j) , in which a first-order transition occurs in 
presence of a broken continuous rotation symmetry. In both cases the same type of 
symmetry-breaking is also known to be possible with a second-order transition. This 
occurs for the standard nearest-neighbour 3-state two-dimensional Potts model, or 
for the standard three-dimensional classical Heisenberg or XY models respectively. 

In fact, the model with many invisible states has a first-order transition for the 



*' E-mail: A. CD. van. Enter®. rug.nl 
*-* E-mail: dajegiu@ginail.com 



typeset using ■pTpT^.ds (Ver.0.9) 



2 



A. CD. ven Enter, G. lacohelli, S. Taati 



same reason the high-g Potts model has a first-order transition. At the transition 
temperature there is coexistence between high-energy phases and a high-entropy 
phase, and between them "free-energy barriers" exist. Such a coexistence can be 
proven by a form of a Peierls-type free-energy-contour argument. For the standard 
Potts model, by now there exists a variety of such proofs,^'' whether by 
Reflection Positivity and Chessboard Estimates, or by a Pirogov-Sinai argument, 
either within a spin description or in a random-cluster version. Typically those proofs 
can be adapted without too much effort to include the model described above. 

Inside standard Potts "free-energy contours" , as were described above, sites exist 
on the border between the ordered and disordered phases and hence they are neither 
ordered nor disordered themselves. They have neither all neighbours different, nor 
all neighbours equal, and as a consequence, in dimension d they lose a free-energy 
fraction per site of order ^ with respect to the free energy of either a disordered 
site (where free energy is purely entropic) or with respect to the free energy of an 
ordered site (whose free energy is purely energetic). At the transition temperature 
the entropy of a disordered site approximately equals the energy of an ordered site. 

In the model we will consider below, the ordered sites need to be in a visible 
colour, whereas disordered sites can have all neighbours either different or invisible. 
It is enough to consider only the two-dimensional version, but this is not essential, 
and the arguments directly generalize to higher dimensions. As the presence of a 
first-order transition in higher-dimensional Potts models is less surprising, the main 
interest seems to be in two dimensions. 

As a further comment we mention that the term "invisible" is actually a bit of a 
misnomer, as at high temperatures the density of "invisible" colours is higher than 
those of the "visible" ones when r > q. Thus most of the colours which appear would 
be the "invisible" ones. 

§2. Main result 

At each site there is a discrete-valued spin variable which can take one out of 
q + r colours, q of which "visible" and r "invisible" . The {q, r)-model then is defined 
by the following (formal) Hamiltonian 

<? 

iJ = - J ^ 6{ai, aj) ^ 6{ai, a) 

{id) a=i 

Here the pairs of nearest-neighbour sites {i,j) live on a lattice of dimension at least 
two. Now we have the following result. 

Theorem 1. For q+r large enough the above model undergoes a first-order transition 
in temperature. At the transition temperature q ordered extremal Gibbs states coexist 
with a disordered extremal Gibbs state. 

Proofs. There are various ways in which one may adapt existing proofs. For example, 
the proof originally due to Kotecky and Shlosman,'^^ later also treated in Refs.fTT]). 
[12]), could be adapted by observing that 
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• Our model has a C-potential (in Georgii's terms), so reflection positivity holds. 

• An ordered bond now will be a bond whose two sites have the same visible 
colour. 

• The "restricted ensemble" for the disordered phase is formed by all configura- 
tions having disordered bonds only, which has an approximate entropy density 
ln((jf + r), when g + r is large. 

Then the arguments used in section 19.3 of Ref. or in Ref. I19|) . using chessboard 
estimates to provide a contour estimate, apply. 

Here we will sketch in some more detail an alternative proof based on the Fortuin- 
Kasteleyn random-cluster description, first derived in Ref. ITTj) . and later treated e.g. 
in Ref. [T3|) . This has the advantage that it extends to values of q and r which need 
not correspond to a spin model interpretation, e.g. q = 1, or q and/or r non-integer. 
To do this we will adapt the ordinary random-cluster representation to include the 
model with invisible states. 

In analogy with the standard Potts model,'^''^ it is possible to rewrite the 
partition function for the (q, r)-Potts model in terms of the partition function for 
a variant of the random-cluster model, which we will call the "r-biased" random 
cluster model. Just as the standard random-cluster model, the r-biased model is a 
correlated bond-percolation model. 

Let G = {S, B) be a finite graph, where S denotes the set of sites, and B the 
set of bonds in the graph. The r-biased random- cluster model on G is given by a 
probability distribution on the sets X C. B. The distribution has three parameters 
0<p<l, (/>0 and r > and is defined by 

in which ko{S,X) denotes the number of isolated vertices of the graph {S,X), 
Ki{S,X) the number of non-singleton connected components of {S,X) and Z^^^(G) 
the partition function. Notice that for r = 0, the model reduces to the standard 
random-cluster model, in which both singleton and non-singleton connected compo- 
nents have weight q. For r > 0, the above model induces a bias towards singleton 
connected components. Namely, the singleton connected components have weight 
{q + r) whereas the non-singleton connected components have weight q. 

Let us now see how the {q, r)-Potts model is related to the r-biased random- 
cluster model. Let 17 be the set of {q, r)-Potts configurations on G. The partition 
function of this model can be rewritten as 

cr&Q 

= ^ Jl e/95(<7,=a,<g) 
<y£n {i,jr}GB 

= E n [^ + 5{a, = a,<q){e^ -I) 

a&a {i,j}GB 
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cr£Q XCB {i,j}eX 
a£QXCB 

where 

7r(a,X) = e'^l^l n [s{{i,j}eX)S{a,=a,<q){l-e-'') + S{{i,j}^X)e-^' . 

The expression above describes a couphng of the {q, r)-Potts distribution on = 
{1, ...q + r}^ and a probabihty distribution on the space {0,1}^ (compare Ref-IS)). 
The marginal of this coupling on the space {0, 1}^ is simply the r-biased random- 
cluster distribution tpp^^q^r with = 1 — e~^. In particular, the weight 7r(cj, X) can 
also be expressed as 

^(a,X)=e^l^l-lp',(a,X). [] [p^ S{{i,j} e X) + {1 - p^) 6{{i,j} ^ X)] , 
where 

Fr = {{(T,X) :ai = aj<q for all {i,j} G X} . 

This expression gives us the model with free boundary conditions. The wired 
boundary conditions in the (g,r)-model correspond to all spins on the boundary 
having the same visible colour. 

The Peierls (free-energy) contour estimate is almost unchanged in comparison 
with the standard case. The only difference is that, (except for the degeneracy term 
of the free contours) the contours of Ref. I17|) satisfy the same Peierls estimate with q 
replaced by g + r. More precisely, in Ref. [13]), in equation (7.59), the adapted Peierls 
estimate for wired contours is of the form 

'I'wilw) < {q + r)-^^ e^l''^"ll , 

whereas the corresponding estimate for free contours, below equation (7.61) of Ref. [T3|) . 
takes the form 

[|7/ll 

^filf) <q{q + r) e^ll'''-^H . 

Convergence of the Peierls estimate will hold once q + r \s large enough, as we then 
have (both for free and wired contours) an estimate of the form 

^(7) < e-^il^ll , 

for a sufficiently large r. 

For a more detailed treatment we refer to Refs. [6]). ll4p . □ 
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§3. Comments and conclusions 



In this note we liave sliown liow the Potts model with many invisible states, 
introduced by Tamura, Tanaka and Kawashima, can be proven to have a first-order 
phase transition, similarly as occurs for the standard high-g' Potts models. The 
transition temperature is asymptotically given by /? « ln{q+r). The proofs, as usual, 
apply for quite high values of q and/or r, the numerical approach of Refs. [20|) . l2ip 
might give a better indication of the values at which the first-order transition first 
occurs. 

We conjecture that the dynamical properties of the Potts model with r invisible 
states, which were considered in Ref. I22|) . have a similar corresponding behaviour 
as occurs in the ordinary Potts model, for which they were rigorously analysed in 
Ref. ED . 
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